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Abstract 

We compute one-loop corrections to the mass matrix of the neutral Higgs bosons of 
the Minimal Supersymmetric Standard Model with explicit CP violation. We use the 
effective potential method, allowing for arbitrary splitting between squark masses. We 
include terms 0(g 2 h 2 ), where g and h stand for electroweak gauge and Yukawa cou- 
plings, respectively. Leading two-loop corrections are taken into account by means of 
appropriately defined running quark masses. 
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1) Introduction 



Supersymmetry is the currently best motivated extension of the Standard Model (SM) of 
particle physics, since it allows to stabilize the gauge hierarchy without getting into conflict 
with electroweak precision data. Among all possible supersymmetric theories, the Minimal 
Supersymmetric Standard Model (MSSM) occupies a special position. It is not only the 
simplest, i.e. most economical, potentially realistic supersymmetric field theory, it also has 
just the right particle content to allow for the unification of all gauge interactions |IJ. Within 
the MSSM, the Higgs sector can be singled out: among all the as yet undetected new particles 
in the MSSM spectrum, the lightest neutral Higgs boson hi is the only one for which a strict 
upper bound on the mass can be given 0, ^ 130 GeV. A good understanding of the Higgs 
sector is therefore of crucial importance for experimental tests of the MSSM. 

It has recently been realized H that explicit CP violation in the mass matrices of third 
generation squarks can induce sizable CP violation in the MSSM Higgs sector through loop 
corrections. Note that CP violating phases for third generation sfermions can be quite large, 
since they contribute to the electric dipole moments of the electron and neutron only at the 
two-loop level f§]. Although a one-loop effect, the induced CP violation in the MSSM Higgs 
sector can be large enough to significantly affect Higgs phenomenology at present 0, |5| and 
future 0, |6| colliders. An accurate treatment of CP violating loop corrections to the MSSM 
Higgs sector is therefore of some importance. 

The first calculations || used diagrammatic methods, and diagonalized the resulting mass 
matrix only approximately. More recently, the effective potential has been used |5[ |7[] to 
compute the 3x3 mass matrix of the neutral Higgs bosons of the MSSM. However, the results 
of ref.|| are not valid for large mass splitting between squark mass eigenstates, while ref. 
does not include contributions from the bottom-sbottom sector, which can be important for 
large ratio of vacuum expectation values tan/5. Here we present a calculation based on the 
full one-loop effective potential, valid for all values of the relevant soft breaking parameters. 
This extends older results [§] where CP was assumed to be conserved. Moreover, unlike 
refs.[|5], [7|, H we include terms 0(g 2 h 2 ), where g and h are a weak gauge and third generation 
Yukawa coupling, respectively. These new terms do not change the spectrum very much, 
but alter CP violating mixing angles by typically 20%. Finally, we absorb leading two-loop 
corrections into appropriately defined running quark masses [|. 

The remainder of this article is organized as follows. In Sec. 2 we present analytical results 
for the 3x3 mass matrix of the neutral Higgs bosons of the MSSM. In the appropriate limit 
we find complete agreement with ref. || , but there are some discrepancies between our results 
and those of ref. J7|. In Sec. 3 we show some numerical results, and compare them with results 
obtained using the formalism of ref.||. Finally, Sec. 4 is devoted to a brief summary and 
conclusions. 



2) Analytical results 

The MSSM contains two Higgs doublets Hi, H 2 , with hypercharges Y(H\) = —Y(H 2 ) = —1/2. 
Here we are only interested in the neutral components, which we write as 

H? = -^(<f> 1 + ia 1 ); H° = ^=(<f) 2 + ta 2 ), (1) 
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where 0i,2 and a\p are real fields. The constant phase £ can be set to zero at tree level, but 
will in general become non-zero once loop corrections are included. 

The mass matrix of the neutral Higgs bosons can be computed from the effective potential 
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where we have allowed the soft breaking parameter m 
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[2 



\rn 



12 



e* 6 * 12 to be complex. We have 



V 



9 2 + a\ - (pi - a\] 



9 



g 2 + g' 2 



(3) 



where the symbols g and g' stand for the SU (2) and U(l)y gauge couplings, respectively. Q in 
eq.(0) is the renormalization scale; the parameters of the tree-level potential, in particular the 
mass parameters mf, m 2 and m\ 2 , are running parameters, taken at scale Q$. The potential 
(0) is then independent of Q , up to two-loop corrections. 

M. is the field-dependent mass matrix of all modes that couple to the Higgs bosons. The 
by far dominant contributions come from third generation quarks and squarks. The (real) 
masses of the former are given by 
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where and h t are the bottom and top Yukawa couplings. The corresponding squark mass 
matrices can be written as 
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Here, and H 2 are given by eqs. (|TJ) while m 2 and m 2 are as in eqs.(|]) and T> has been 
defined in eqs.(§). In eqs.(|5|) m~, m~ and m~ are real soft breaking parameters, Aj, and A t 

are complex soft breaking parameters, and fi is the complex supersymmetric Higgs(ino) mass 
parameter. The eigenvalues of the mass matrices (|5|) are: 
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The calculation proceeds by plugging the field-dependent eigenvalues (f|) and (§) into the 
potential (0); here the complex scalar squarks and Dirac fermion quarks enter with overall 
factors +2 and —4, respectively. The mass matrix of the Higgs bosons (at vanishing external 
momentum) is given by the matrix of second derivatives of this potential, computed at its 
minimum. In order to make sure that we are indeed in the minimum of the potential, we solve 
the stationarity relations, i.e. set the first derivatives of the potential to zero. This allows us 
to, e.g., express m 2 , m\ and m 2 2 sin(£ + 6 12 ) as functions of the vacuum expectation values 
(vevs) and the remaining parameters appearing in the loop-corrected Higgs potential. Note 
that the equations dVmggs/dai = and dVmggs/ 'da 2 = are linearly dependent, i.e. lead to 
only one constraint on parameters, if we demand that (ai) = (02) = 0; the remaining vevs are 
defined through 
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This leads to the following expression for the re-phasing invariant sum £ + # 12 : 
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where 



In eq.( 
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we have introduced the quantities 
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which describe the amount of CP violation in the squark mass matrices. Note that A^ remains 

since this implies ^sm(A t /j,e 1 ^) — > 0. Most tree-level analyses use the 
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convention £ = 0; in this case we could set it to zero in the right-hand side of eq.(|§|), to one- 
loop order. We kept it in order to illustrate that only the phases of the re-phasing invariant 
quantities A t fie^ and A&//e^ have physical meaning. 

The mass matrix of the neutral Higgs bosons can now be computed from the matrix 
of second derivatives of the potential (§), where (after taking the derivatives) m\, m\ and 
m 



12 



sin(£ + 6*12) are determined by the stationarity conditions. We find that the state G 



a± cos P — a.2 sin f3 is massless; it describes the would-be Goldstone mode that gets "eaten" by 
the longitudinal Z boson. We are thus left with a squared mass matrix M. 2 H for the three 
states a = disin/3 + a 2 cos/?, and 02- This matrix is real and symmetric, i.e. it has 6 
independent entries. The diagonal entry for a reads: 
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where A^ and A^ are as in eq.([H]), and we have introduced the function 



The quantity m\ in eq.flTT]) is given by 
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where 
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the function / has been defined in eq.@. 

If we consider m\, the values of the soft breaking parameters and \x to be inputs, eqs.([E 
and (O) can be combined to give an explicit expression for the induced phase £ of the vev: 
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It should be emphasized that £ and 9\2 are not separately physical quantities; only their sum is 
re-phasing invariant. Nevertheless eq.(|i~5D can be useful. Note that the leading Q dependence 
in the numerator cancels between the explicit terms oc logQo contained in F(m~^,m~^) and 
F(m^, m| ), and the RG-induced running of 3rn(mf 2 )- If one sets the phase of the Higgs field 
at some scale, eq.fllS]) shows that it will remain scale-independent, at least to 1-loop order. 
In particular, £ can be set to zero, as pointed out by Pilaftsis ||. Eq . (|i"5"D then determines 
the loop-induced phase of m\ 2 (called a counter-term in ref.[H]). On the other hand, it is also 
possible to set Qyi = 0. The price one has to pay is that the phase of the Higgs fields will show 
a strong scale dependence even after explicit 1-loop corrections to the Higgs potential have 
been added. 

We will see below that most of the Q dependence of the loop corrections to the Higgs mass 
matrix can be absorbed into the parameter m^. Within m\ itself, the Qq dependence largely 
cancels between the running of m\ 2 and the explicit logQo dependence of f(m 2 ), eq.(|), just 
as in the CP conserving case ||. Note that our itla differs from M a defined in ref.||, since 
our rriA only includes corrections that are nonzero in the limit of exact CP invariance, whereas 
M% corresponds to our A^^| aa of eq. (|TT|) . In the remainder of this paper we will consider m\ 
as well as the re-phasing invariant quantities A t fie^ and A^fie^ to be input parameters. The 
results are then independent of the convention adopted for # 12 . 



The CP violating entries of the mass matrix, which mix a with 0i and 02, are: 
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As noted earlier, the size of these entries is controlled by A t ~ and A b . In addition we have 
introduced the dimensionless quantities 
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As in eqs.fllOD, phases only appear in the re-phasing invariant combinations Atfie^ and A b fie 1 ^. 
The terms proportional to g 2 [defined in eqs.(§)], X t or X b are mixed gauge- Yukawa contri- 
butions, which were neglected in refs.[]5], |7J; note that no corrections of this kind appear in 
eqs.(||) and (pi]). We do not include pure gauge, 0(g 4 ) corrections, since there are many 
additional corrections of this order from first and second generation sfermions as well as from 
loops involving gauge and Higgs bosons and their superpartners. 

While eqs. (|T6|) are pure loop corrections, the remaining entries of the Higgs boson mass 
matrix also receive tree-level contributions: 
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As noted earlier, the explicit dependence on the renormalization scale Qo has mostly been 
absorbed into m\. The only exceptions are terms proportional to the combination g 2 of 
electroweak gauge couplings, eqs.(0). They come from the wave function renormalization of the 
Higgs fields, which leads to a logarithmic Qq dependence of the tree-level contributions g 2 ((j>i) 2 
( in MW^), -^ 2 (0i)(0 2 ) (in M 2 H \ M2 ) and <? 2 (0 2 ) 2 (in M 2 H \^), respectively. Following 
ref. || we define these vevs at scale Qq = m t , i.e. we use Qq = m t in the explicitly Q dependent 
terms in eqs. (|l8D . 

In order to compare our results with those of Demir 0, we have to set g — g' — hb — 
in the loop corrections. Our results then become similar to his if we identify our m\ with 
his M\. However, his CP-odd quantity A seems to be too small by a factor of 2, and the 
coefficient of g(m 2 ^,m 2 ^) in -M-hI^^ differs if CP is violated. 

In order to compare our results with those of Pilaftsis and Wagner || , we have to set most 
corrections involving gauge interactions to zero, the only exception being the terms oc g 2 log Qq. 
In addition, we have to take the limit of small squark mass splitting, which implies 
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where M| USY = + m | 2 ) After these replacements we find complete agreement with 
the results of ref.[|], except for *M^| aa ; however, this difference can be absorbed into the input 
parameter m 2 A) i.e. it has no physical significance. 

Pilaftsis and Wagner also include logarithmically enhanced two-loop corrections 0(g 2 h 4 ) 
and 0(h 6 ), where g s is the strong gauge coupling. As shown in refs.|| leading QCD 
corrections can be included by interpreting the quark masses and Yukawa couplings appearing 
in eqs.(|ll|), fll3D , (|TJ) and (|18D to be running parameters. For example, 
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where we have assumed Q < Msusy, i-e. only standard QCD contributes; a pole mass 
777,P ole = 173 GeV corresponds to rnt{m t ) = 165 GeV. The corresponding expression for the 
running Yukawa coupling follows from eq.(|4]) and the observation that, to one-loop order, QCD 
corrections do not contribute to the running of (fa)- Corrections 0(h 6 ) to the Higgs mass ma- 
trix come from the one-loop running of the Yukawa couplings, as well as (for Q < Msusy) from 
the two-loop running of the quartic Higgs coupling(s) in the effective non-supersymmetric the- 
ory. Finally, as pointed out in ref. 0, one can absorb significant 0(g 2 h A ) corrections that are 
not enhanced by large logarithms by including gluino-stop loop corrections to the running 
top mass at scale Q ~ Msusy- In order to (approximately) include all these corrections, we 
introduce three different top masses that appear in various contributions to Ai 2 H : 

9 9 

• The leading-log term oc \h t \ m t log H 4 * 2 should be computed with a running top mass 
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at the intermediate scale Qf nt = m t MsvsY- 
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In terms proportional to powers of the CP-odd quantity Af, one should use the top mass 
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Finally, in the remaining terms the top mass should be interpreted as 
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The quantity m t ^igh appearing in eqs.( [22"| , [2"3"|) is given by 
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Here 6i and fa are the angles needed to diagonalize the stop mass matrix (|5a|), i.e. t± = 
cos^f ti + sin^e - ^* tR. m g is the gluino mass, which can be complex, and B is the Passarino- 
Veltman two-point function. Eq.([24|) generalizes corresponding results of ref. to the case 
with CP-odd phases; note that the quantity e'^'mj is re-phasing invariant. Eqs.(|2l|)-(24) 
can be extended straightforwardly to the bottom-sbottom sector.fj] 



3) Numerical examples 

We are now ready to present some numerical results. It is clear from eqs.flTED and ([It]) that 
loop-induced CP violation in the Higgs sector can only be large if both and \A t \ (or \Ab\, 
if tan/5 3> 1) are sizable @. We therefore choose \A t \ = \A b \ = = 2m^. For definiteness we 
only present results for fixed tan/3 = 10, and real and positive gluino mass. The value of tan/3 
is not very important, unless it reaches ~ m t jm^ where b — b contributions can be important. 
Since for the given moderate value of tan/5 contributions from the (s)bottom sector are still 
quite small, our result are not sensitive to m-^ and we therefore fix = and also take 
equal phases for A t and A b . However, we allow different values for the soft breaking masses 
of SU(2) doublet and singlet squarks, m^ ^ m^j. This allows us to study scenarios with large 

t\ — t2 mass difference even if 3> m t . Note that renormalization group effects do in fact 
produce significant differences between and even if they are equal at some very large 
energy scale |p] . Finally, we take m t (m t ) = 165 GeV and mb(mb) = 4.2 GeV. 



We first study a scenario with = m^. Here we expect good agreement with results of 
ref. || as long as 3> m t , if we "switch off" the new contributions oc g 2 h 2 presented in Sec. 2 
as well as the SUSY-QCD correction (|24]). This expectation is borne out by Fig. 1, which 
shows the mass of the lightest neutral Higgs boson (left panel) as well as the CP violating Higgs 

^In eq.(|24|) we have assumed that mt < nig + mj . If this is not the case, the loop functions will develop 
imaginary parts. To one-loop order, only the real part of the correction contributes. 
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mixing angle 02 (right panel) as a function of mg, for itla = 1 TeV and 8 c s = arg(A^e l * 
7r/2. Defining the orthogonal Higgs mixing matrix H through hi = 0^<pj, where hi dene 
the physical Higgs bosons (mass eigenstates) and we have set 03 = a, the ai are given by 



mm 



oKr+\o. 



12 I 



o{l\ 2 + \o%\ 2 



\o H \ 

1^3 I 



(25) 



If CP is conserved, hi is either purely CP-odd (O^ = 1), or purely CP-even ( Of{ 



+ 



1), in which case a,- = 0. 




Figure 1: The mass of the lightest neutral Higgs boson (left) and the CP violating mixing angle of 
the second-lightest neutral Higgs boson (right) as a function of the soft breaking mass m~ of SU(2) 
doublet squarks. The values of the other free parameters are: m~ = m~ = m^, \A t \ = \Ab\ = |/i| = 
2m~, tan/3 = 10, tua = 1 TeV, m~„ = 1 TeV and squark phase 6 e s = tt/2. The dotted lines have 
been obtained using the expressions of ref.||, while the dashed and solid curves show our results 
without and with 0(g 2 h 2 ) and 2-loop SUSY-QCD corrections, respectively. 

The dotted curves in Fig. 1 have been obtained using the expressions of ref. . As expected, 
for large we find very good agreement with the dashed curves, which show our results with 
the new 0(g 2 h 2 ) and 2-loop SUSY-QCD contributions switched off. For this comparison we 
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have chosen the parameter of ref.0 to agree with the (a, a) element (|11|) of the Higgs 
mass matrix. We then find that the expressions of ref.[f| reproduce the mass of the lightest 
neutral Higgs boson very well even for quite small values of mg. This is somewhat surprising, 
since for the smallest experimentally allowed value of mx, defined through the requirement 



m t x > 90 GeV fL2| , the stop mass splitting is very large, m ( - 2 ~ 477^. On the other hand, the 
two predictions for a 2 do start to deviate significantly for m-x ^ 500 GeV, which corresponds 

to mi 2 ^2m^. Note that «i is very small in this example, 0(1O~ 3 ), due to the large difference 
between and m hl . 

The solid curves show our full results, including all contributions described in Sec. 2; we 
choose a gluino mass of 1 TeV. The new contributions to are not very large, although 
the reduction by 4 GeV for the smallest allowed value of is not entirely negligible. In 
contrast, these new contributions reduce a 2 by ~ 20% for all values of m^. Contributions of 
this magnitude should certainly be taken into account when translating observed (bounds on) 
CP violating effects into (bounds on the) values of the fundamental soft breaking parameters. 




Figure 2: The mass of the lightest neutral Higgs boson (left) and the CP violating mixing angle of 
this Higgs boson (right) as a function of m~. Notation and parameters are as in Fig. 1, except that 
we have taken m~ = m~ = '^ m Q an d rriA = 100 GeV. 

In Fig. 2 we show results for a scenario with = 2rrig, so that m~ t2 > 1vn~ tl for all values 
of m-x. We nevertheless find good agreement between our calculation of and the results 
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of ref. ||. The reason is that we now have chosen uia = 100 GeV, in the range of interest to 
current LEP experiments. In the absence of CP violation we would then have ~ tua H, 
since tan 2 (3 3> 1 in our example. In the presence of CP violation in the squark sector the 
corrections (|ll]) to the (a, a) element of the Higgs mass matrix are negative; is reduced 
further by mixing between CP-even and CP-odd Higgs bosons. However, these effects remain 
quite small as far as the Higgs spectrum is concerned, and therefore need not be computed 
with very high accuracy. 

In contrast, the three predictions of a± now differ significantly over the entire allowed 
range of m^. In the given case we observe a partial cancellation between effects due to squark 
mass splitting and the new 0(g 2 h 2 ) and 2-loop SUSY-QCD contributions. As a result the 
prediction from ref. |J agrees more closely with our full result than with the approximation 
that only includes the interactions also included in ref. []5| . We see that the new mixed gauge- 
Yukawa and SUSY-QCD corrections again reduce CP violation in the Higgs sector by 10 to 
20%. Note, however, that this loop-induced CP violation in the Higgs sector can still be very 
large, i.e. a\ ~ 0(1) is possible here. 

Finally, in Fig. 3 we show that oti can in fact reach its theoretical maximum of unity. In 
this figure we keep all dimensionful parameters fixed, and show m/ ll and a\ as a function of 
the phase 6 e g appearing in the squark mass matrices. The dependence of m/ ll on this phase is 
mild. Note that the lightest Higgs mass now takes its maximal value for maximal CP violation 
in the squark sector, 6 e s = n/2. In contrast, the CP violating Higgs mixing angle a\ becomes 
maximal for intermediate values, 9 e g ~ 0.9 and 8 e s ~ 2.2. In between these two values, hi is 
mostly CP-odd, i.e. > l/y/2; in the remaining region hi is dominated by its CP-even 

components. 

Note that a\ is actually quite small, ~ 0.1, for "maximal" CP violation in the squark 
sector. The reason is that the (a, <pi) element (|16a|) of the Higgs mass matrix is small here. 
The leading contribution to this element is proportional to the quantity R t of eq. (|17b| ), which 
is small if A t fie % ^ is purely imaginary, since the contribution oc \fi\ 2 is suppressed for our 
choice tan/? = 10. As a result, this element of the Higgs mass matrix goes through zero for 
6 c q ~ 0.537T. The quantity R' t remains large in this region of parameter space, leading to a 
sizable (a, 02) element of the Higgs mass matrix. However, its contribution to «i is suppressed 
by the relatively large difference between the diagonal (a, a) and (02,02) elements. 



4) Summary and conclusions 

In this paper we calculated quark and squark loop corrections to the mass matrix of the 
neutral Higgs bosons of the MSSM using the effective potential method. We allowed for CP 
violating phases in the squark sector as well as arbitrary splitting between squark masses. We 
also for the first time included mixed weak gauge- Yukawa contributions, as well as leading 
2-loop SUSY-QCD contributions, in the calculation of induced CP violation in the Higgs 
sector; the latter have been absorbed into properly defined running quark masses, following 
ref. |2j . Our formulas can be used to accurately calculate the masses and mixing angles of the 
MSSM Higgs bosons for all values of the (possibly complex) parameters describing the third 
generation squark mass matrices. 

When comparing our results with earlier expressions || that are valid for relatively small 
squark mass splitting we found surprisingly good agreement for the Higgs spectra even in cases 
with very large squark mass splitting. However, we also found ~ 30% discrepancies in the 
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m^=m^=m^=0.5 TeV, m A =0.1 TeV 




Figure 3: The mass of the lightest neutral Higgs boson (dashed, referring to the scale at the left) 
and the CP violating mixing angle of this Higgs boson (solid, referring to the scale at the right) as a 
function of the effective squark phase 6 e g = arg(^4 t/ ue^). We have fixed m~ = m~ = m~ = m~ g = 500 
GeV and txia = 100 GeV; the values of the other parameters are as in Fig. 1. 

predictions of CP violating Higgs mixing angles if the difference between squark masses is large. 
Similarly, we found that the mixed weak gauge- Yukawa and 2-loop SUSY-QCD contributions 
do not change the predictions for the Higgs spectrum very much, but can change the prediction 
for CP violating Higgs mixing angles by ~ 20%. 

Finally, we showed that maximal mixing between CP-even and CP-odd Higgs states re- 
mains possible for reasonable choices of the free parameters. Note that this usually does not 
happen when CP violation in the squark sector is maximal, i.e. CP violating Higgs mixing 
angles reach their maximum for intermediate values of the effective CP violating phase in the 
scalar top mass matrix (away from 7r/2). Since the loop-induced CP violation in the MSSM 
Higgs sector can be very large, an accurate treatment of this effect, as described in this paper, 
is required for the interpretation of searches for Higgs bosons at LEP and elsewhere. 
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